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ABSTRACT 

The ACS Survey of Galactic Globular Clusters is a Hubble Space Telescope 
{HST) Treasury program designed to provide a new large, deep and homo- 
geneous photometric database. Based on observations from this program, we 
have measured precise relative ages for a sample of 64 Galactic globular clusters 
by comparing the relative position of the clusters' main sequence turn offs, us- 
ing main— sequence fitting to cross— compare clusters within the sample. This 
method provides relative ages to a formal precision of 2-7%. We demonstrate 
that the calculated relative ages are independent of the choice of theoretical 
model. Wc find that the Galactic globular cluster sample can be divided into 
two groups - a population of old clusters with an age dispersion of ~5% and no 
age— metallicity relation, and a group of younger clusters with an age— metallicity 
relation similar to that of the globular clusters associated with the Sagittarius 
dwarf galaxy. 

These results are consistent with the Milky Way halo having formed in two 
phases or processes. The first one would be compatible with a rapid (<0.8 Gyr) 
assembling process of the halo, in which the clusters in the old group were formed. 
Also these clusters could have been formed before reionization in dwarf galaxies 
that would later merge to build the Milky Way halo as predicted by ACDM cos- 
mology. However, the galactocentric metallicity gradient shown by these clusters 
seems difficult to reconcile with the latter. As for the younger clusters, it is very 
tempting to argue that their origin is related to their formation within Milky Way 
satellite galaxies that were later accreted, but the origin of the age— metaUicity 
relation remains unclear. 
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Subject headings: Galaxy: evolution — Galaxy: formation — globular clusters: 
general 



Introduction 



One of the keys to understanding the structure and evolution of the Milky Way is the 
ability to div ide its stars and clusters into separat e Gala ctic Populations as first hinted at by 
Port I (Il926l ) and formally introduced by iBaade I (Il944l ). In the early 1950s, Sandage pub- 
lished the first deep color— m agnitude diagrams (CMDs) of the northern G alactic globular 
clusters (GGCs) M3 and M92 (ISandage lll953l : lArp. Baum. fc Sandage Ill953l ) and discovered 
the faint sequence of stars similar to the main sequence of Population I. It was almost im- 
mediately realized that GGCs are ancient compared to the stars in the Solar Neighborhood. 

Over the last half century, it has become clear that GGCs are astronomical fossils 
providing important information on the formation and evolution of the Milky Way. One 
characteristic of GGCs that makes them so valuable is that they are the oldest Galactic 
objects for which reliable ages can be measured. 

GGC ages are of great interest to a variety of cosmological and comogonical issues 
since their absolute ages place a lower limit on the age of the Universe. Unfortunately, 
absolute age measurements are still compromised by a number of uncertainties, particularly 
GGC distances and foreground reddenings as well as metallicities, bolometric corrections 
and some poorly known aspects of the stellar evolution. Relative ages are almost (although 
not totally) free of these effects and can still reveal fundamental information about Galaxy 
formation mechanisms and time scales. For this reason they have become of paramount 
importance in the context of the study of GGC. 

The study of GGC ages has evolved apace with improvements in instrumentation, in 
our theoretical understanding of stellar evolution and with the advance of computational 
facilities. Reliable age determination requires deep photometry that reaches at least the 
main— sequence turn off (MSTO), together with a good understanding of the stellar evolution 
theoretical modeling. During the first half of the 20*^ century, the available instrumentation 



^ Based on observations with the NASA/ESA Hubble Space Telescope, obtained at the Space Telescope 
Science Institute, which is operated by AURA, Inc., under NASA contract NAS 5-26555, under program 
GO-10775 (PI: Sarajedini). 
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did not allow observations deeper than the horizontal branch of most GGCs. Over the past 
three decades, however, low quantum efficiency non— linear photographic imaging has given 
way to digital imaging using high quantum efficiency linear CCDs. This has resulted in ever 
deeper and more precise cluster CMDs, that reached the MSTO and allow ed the app l icatio n 



of age da t ing te c hniques. First relat i ve ag e stu dies were carried out bv iGratton I (119851 ) 



Peterson I (119871 ) , ISarajedini fc King I (119891 ) and ISarajedini fc Demaraue I (119901). as 



sentative examples. However, uncertainties in both the input physics ( IChaboyer et al. Ill998l ) 



repre 



and the color— Tpff tra nsformation of the theoretical models into the observational plane 



(iBuonanno et al. Ill998l ) constitute a significant source of errors in the final relative age deter- 



mination, particularly when ages of clusters with different metallicities are compared. In the 



^^„v. , , „ „ w ^. rt^-^, w . ^ 

with each technological or theoretical advance (e. s. Chabover. Demaraue. & Sarai 


edini 


19961: Richer et a 


. 199d: Stetson et al. 19961: Salaris & Weiss 


1998: Buonanno et al. 


1998; 


Rosenberg et al. 


I999I: 


VandenBerg 


2OO0I: De Angeh et al. 2005h. 



Relative ages can be estimated to high precision by measuring the position of the 



MSTO relative to CMP fea tures that have httle or no age dependence (IStetson et al 



19961 : ISarajedini et al. 1 119971 ). However, a key aspect of any relative age experiment is 
having a large homogenous database of globular cluster photometry. Such databases have 
only recently appeared in the literature but have significantly advanced our understand- 

,nd evolution. The largest dedicated eff orts are those of 



Rosenberg et al. 


( 


2000a 


b), 


1999 


), and 


De Angeli et al. 



(l2000alJbl') . a grou nd-based survey of 35 nearby G GCs (IRosenberg et al. 



(I2OO5I ) , an analysis of 55 clusters from iPiotto et al. I (120021 ) 's 



HST snapshot catalog. 



These and other recent relative age studies (IBuonanno et al. I Il998l : ISalaris fc Weiss 



19981 : [Rosenberg et al. I Il999l : IVandenBerg I l2000l : iDe Angeli et al. I l2005l ) consistently point 
towards a scenario in which metal— poor clusters are largely coeval but metal— rich clusters 
show some age dispersion, ~1.0 Gyr (rms), and a possible age— metallicity relation. No 
correlation between age and Galactocentric distance has been conclusively demonstrated. 

In this paper, we present an analys is of a uniform se t of HS T/ACS data for 64 GGC^. 



This dataset, previously introduced by ISarajedini et al. I (120071 ) . has produced well-defined 



CMDs for all of the target clusters from the tip of the re d giant branch (RGB) to at least 
~6.5 magnitudes below the MSTO (lAnderson et al. 1120081 ). The determination of GGC ages 
was a principle goal of the HST/ACS survey and exposures times were defined to maximize 



iSaraiedini et al. I (|2007r ) presented observations of 65 GCs, but Pal 2 has been left out of this study due 
to its high differential reddening 
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the S/N ratio at the MSTO. This makes the survey an excellent resource for the evaluation 
of relative ages. 



Due to the difficulty in determining the position of the MSTQj to the high precision 
required for age studies, many different m ethods have been developed in previous relative 
ages works (see iMeissner. fc Weiss I l2006l . for a complete discussion). The two most used 
relative age indicators are those based on the color difference between the MSTO and the 
RGB at a given magnitude level (horizontal method), and the magnitude difference b e tween 
the MSTO and the zero— age HEg (vertical method). IVandenBerg. Bolte. fc Stetson I (llQQOl ) 
were pioneers in measuring GGC relative ages making use of the horizontal method. They 
essentially derived relative ages comparing the color difference between the MSTO and the 
location of a well defined point in the lower RGB, located 2.5 magnitudes brighter than a 
point in the upper main— sequence (MS) that is 0.05 magnitudes redder than the MSTO. 
However, horizontal method results are sensitive to the assumed value of the mixing length 
parameter and to the color— Tog transformation of the theoretical models into the observa- 
tional plane. As a result, relative ages determined using the horizontal method are model 
dependent. Potentially, more reliable age indicators are those related to the brightness of 
the MSTO, in particular the vertica l method. First surve ys of GGC relat ive ages using the 
vertical method were carried out by iGratton I (119851 ) and iPeterson I (119871 ). They measured 
relative ages from a sample of 26 and 4 1 GGC s, respectively, collected from the available 
literature. Later on, ISarajedini fc King I (119891 ) used the brightness difference between the 
HB and the MSTO to estimate the ages of 31 GGCs. Nevertheless, accurate determination 
of the zero— age HB is delicate, in particular for the clusters at the extreme boundaries of 
the metallicity distribution, which generally have a red (or blue) HB and no RR Lyrae stars. 
Moreover, the dependence of the zero— age HB luminosity on the cluster metallicity is still 
controversial. 

In this paper, we use an alternative approach. Our data reaches at least ~6.5 magnitudes 
below the MSTO, providing a very well defined MS. This allows, for the first time, the use of 
relative MS fitting for a large number of clusters. As it will be shown later, the relative MS 
fitting procedure used in this study does essentially the following. The MS of two clusters, 
with similar metallicities, are superimposed in the CMD by shifting one of the cluster's 
mean ridge line in both color and magnitude (explained in detail in section H]). This way. 



The MSTO point is the current position in the CMD where stars in the cluster start burning hydro- 
gen in a shell and are leaving the main— sequence. The empirical MSTO is then the bluest point on the 
main— sequence— to— subgiant— branch portion of the cluster's CMD. 

The zero— age HB is the position in the CMD where stars in the cluster start burning helium into carbon 
and oxygen in a helium rich core and hydrogen into helium in a shell. 
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any differences between the two cluster's distance and reddening are compensated, and an 
intrinsic MSTO magnitude difference is obtained. Tfiis relative MSTO brightness is then 
used to derive relative ages. This allows us to use what is, in essence, an improved version 
of the traditional vertical method that substitutes the well-defined MS for the contentious 
zero— age HB. 

We also take the additional step of comparing our observational measures t o updated 
theore tical models. We have primarily used the most recent state— of— art models of lDotter et al. 
( 120071 ) (hereafter D07), transformed to the observational plane with an er apirical Color— T^ff 



transformation. Howe ver, we also cq r apare our results to the r aodels of iPietrinferni et al. 



fl2004[ ) (hereafter P04). lBertem et al. I fll994l ) (hereafter B94) and lGirardi et al. I (120001 ) (here- 



after GOO). We find that the various theoretical models agree closely on the MSTO absolute 
magnitude and color, resulting in negligible model dependence. 

This paper is organized as follows. Section [2] describes the observations and data re- 
duction pipeline, and Section [3] presents the globular cluster database. The mean ridge line 
fitting and MS— fitting procedures used in this study are described in Section HI The analysis 
of the uncertainties is presented in Section O Finally, relative age results and a discussion are 
presented in Sections [6] and Ui respectively. The paper ends with a summary of conclusions 
(Section E]). 



2. Observations and data reduction 

Observations were performed with the ACS/WFC instrument on board the HST. A 
total of 65 GGCs were observed in the F606W (~ V) and F8UW (~ /) fihers. Due to its 
high differential reddening. Pal 2 is not considered here. 



Details on the observations, da ta reduction and photometry are given in lSarajedini et al. 



(120071 ) and [Anderson et al. I (120081 ). For most globular clusters, the observation and data 
reduction pipeline produces ~12 magnitudes of precise photometry, extending from nearly 
the tip of the RGB to several magnitudes below the MSTO. In general, the photometry 
reaches a depth corresponding to an approximate stellar mass of ~ 0.2 m©. 

This uniform and deep photometry offers a database of unprecedented quality, opening a 
window for new data analysis approaches. In particular, cluster MSs are so well defined that 
it is possible, for the first time, to apply relative MS— fitting techniques to a large number 
of clusters within the same database. 
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The database 



The position of the MSTO depends not only upon age but also upon metallicity. This 
effect requires the use of a large and homogeneous metallicity database in order to determine 
accurate relative ages. To account for this, we have adopted the rn e tallici ties listed in the 
Rutledge etaTI Jl997 ) catal ojjl over both the ICarretta fc Gratton I Jl997h (hereafter CG) 
and the IZinn fc West I (119841 ) (hereafter ZW) metallicity scales. This catalog is maybe not 
as precise as a true compilation of metallicities based on high dispe rsion spectra, bu t it ha s 
the the advantage of being homogeneous. For clusters not listed in iRutledge et al. I (Il997l ). 
the original metallicity listed in th e ZW catalog has been transformed to the CG metallicity 
scale using the following equation ( Carretta et al. II2OO1I ): 



[Fe/H]cG = 0.61 + 3.04[Fe/i/]zw + 1.981 [Fe/i/]|w + 0.532[Fe/i/ 



3 

ZW 



Lynga 7, E3, and Pal 1 metallicities are not listed in the previous catalogs. Lynga7 
and E3 metallic i ties in the CG scale hav e been calculated from the ZW metallicities of 
Tavarez &: Friel I (I19951) and [Harris Jl996h . respectively. For Pal 1, we used the measure of 



Rosenberg et al. I (119981 ) 



Table [T] presents the list of our target GGCs, grouped by metallicity. Column 3 gives 
the cluster identity while columns 4 and 5 present the adopted metallicities in the ZW and 
CG sca l es, res pectively. Column 6 lists the Galactocentric distance (roc, in kpc), taken from 



Harris I (Il996f ). 



4. Mean ridge lines and main sequence fitting method 

In this section, we describe the general properties of the cluster mean ridge line fitting 
procedure. We then use the derived ridge lines to perform relative MS— fitting between 
clusters. 



' Rutledge et al. ( 1997h presented a large homogeneous metallicity database (over 71 GGCs) based 
calcium index (Ca II Triplet) measurements. 
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4.1. Mean ridge lines and main sequence turn-off apparent magnitude 



In order to determine accurate relative ages for each cluster, a precise determination of 
the MS and RGB mean ridge line is required. With this aim, we have created a new software 
package for the ACS program. It is worth mentioning that mean ridge lines are determined 
in the observational plane {F606W-F8UW , F8UW) instead of {F60QW-F8UW, FQOQW), 
as the SGB appears more vertical in the CMD and the following procedure produces more 
accurate results. For the rest of the analysis, the {FGOGW-FSliW , F606W) plane is used. 

For each globular cluster, the mean ridge line is determined in two steps (Figured]). 
First, color histograms are computed for stars grouped in 0.4 magnitude {F814W) wide 
bins. The value of 0.4 magnitudes was determined experimentally as a compromise between 
two effects - the low resolution produced by larger bins and the noisy ridge line produced by 
smaller bins. Color histograms are constructed along the RGB and MS with a step of 0.04 
magnitudes (a tenth of the bin size) in F814W. That is, a moving bin of fixed width (0.4 
magnitudes), with successive steps of 0.04 magnitudes, is used. A preliminary mean ridge 
line is then derived with {F606W — F814W) colors set at each histogram's maximum, and 
the F814:W magnitude of the bin center. 

Once the preliminary mean ridge line has been determined, the fitting process is iterated 
three more times. However, we now use rotated histograms. For each F814W value, the 
tangential angle of the mean ridge line is determined. This angle is used to rotate the refer- 
ence system on which the histogram is computed. The X— axis of the new reference system 
is perpendicular to the preliminary mean ridge line at the considered F814W magnitude. 
The histogram is then computed perpendicular to the preliminary mean ridge line and con- 
sidering stars in the new reference system's Y— axis interval [-0.04, 0.04]. The histogram's 
maximum is determined, defining a point of the new mean ridge line in the rotated reference 
system, {max ,0). Note that this point shows the correction to be applied to the preliminary 
mean ridge line. The coordinates of this point are de — rotated and its location in the original 
reference system, the CMD, is calculated. The process is repeated for all the F814:W values 
to derive the revised mean ridge line, which is then smoothed to derive the final mean ridge 
line. Figure [2] shows some examples of mean ridge line fitting results. Four selected clusters 
covering the metallicity range are shown. 



For cl usters with multiple stellar populations, such as NGC 18 51 (IMilone et al. 
NGC 2808 JPiotto et al I boOTh . M54 Jsiegel et al. IboOTh and cuCen Lee et al. Ill999 



2008a') 



Pancino et al. 



2000l : iBedin et al. II2004I ). the mean ridge line fitting procedure fits the dominant population 
(see Figured]). 



The empirical MSTO is the bluest point on the mean ridge line. It is worth mentioning 
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that the final mean ridge fine consists of a discrete number of points, it is not a continuous 
fine. For this reason, in order to determine the MSTO position, a continuous sphne has been 
fitted to the main— sequence— to— subgiant— branch portion of the final mean ridge line. Its 
bluest point is then adopted as the empirical MSTO. As an example, the position of the 
MSTO is marked with an open circle in Figure [2l In order to increase the precision of the 
MSTO position, the mean ridge lines have been determined ten times for each cluster, with 
the starting point of the histograms offset in steps of 0.004 magnitudes. The MSTO color and 
apparent magnitude have been determined as the mean of the ten different measurements. 

This method provides a consistently defined mean ridge line, minimizing the effects of 
the binary population and differential reddening (see section [5]). Additionally, the rotated 
coordinate system ensures that the best fit is derived around the steeply curved MSTO 
and sub giant branch region of the CMD. For high metallicity clusters, in particular, the sub 
giant branch is nearly horizontal and the more commonly used mean ridge line determination 
method (drawing an approximate ridge line, rendering the sequence to a nearly vertical line 
and then taking histograms in small ranges of magnitude) lacks the resolution around the 
sub giant branch needed for precise age measurement. 



4.2. Main sequence fitting method 

The CMD location of the faint MS of a GGC is independent of its age, but highly 
dependent upon its metallicity. For this reason, we have divided the HST/ACS sample 
into six metallicity groups (see Table [1]) using the following divisions: [Fe/H] = -0.3, -0.8, 
-1.1, -1.3, -1.5, -1,8 and -2.3. For each group, a reference cluster has been chosen: NGC 
6304, NGC 6723, NGC 6981, NGC 6681, NGC 6101 and NGC4590, respectively. Reference 
clusters have been selected to have low differential reddening, low field star contamina- 
tion, and a well defined MS. Using the derived mean ridge lines, we performed MS— and 
RGB— fitting between each cluster in each metallicity group and the corresponding reference 
cluster. Figure [3] shows examples of the fitting, in which the mean ridge lines of clus- 
ters with —1.3 < [Fe/H]QQ < —1.1 (sohd lines) have been shifted in both magnitude and 
color to fit the reference cluster, NGC 6981 (dashed line). The fitting was performed in a 
least squares fashion and taking into account two CMD regions that have little dependence 
upon cluster age. The two intervals are [(MpQQg-y^^ — 2.5) < Mpeoew < (^reoew ~ ■^•5)] 
[(^F606W + 1-5) < Afpeoew < (^f606W + 3-0)], Mjeoew being the magnitude of the refer- 
ence cluster's MSTO. These two intervals have been defined making use of D07 theoretical 
isochrones. In particular, D07 isochrones with similar metallicity and different ages were 
superimposed on the same CMD. Based on visual inspection, these two particular regions 
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were found to have little dependence upon cluster age, and were adopted as the optimum 
intervals for the MS fitting procedure. These regions are shaded in Figure [31 This fitting 
procedure produces unequivocal results, with no degeneracy between distance modulus and 
reddening. 

The MS— fitting provides relative MSTO magnitudes for clusters inside each metallicity 
group. In order to calculate relative MSTOs between two consecutive metallicity groups, we 
applied the MS— fitting technique between the corresponding reference clusters. This was 
repeated for every two adjacent metallicity groups, producing a unified sample of relative 
MSTO magnitudes for all the analyzed clusters. 

In order to compare relative MSTO measures to model predictions, we had to calculate 
the absolute magnitudes of the turnoffs. As the relative MSTO magnitudes are known, the 
MSTO absolute magnitude is only needed for one cluster. We chose the well-measured NGC 
6752 cluster, adoptin g a subdwarf-based distance modulus of (m — M)y = 13.24 ± 0.08 from 



Gratton et al. I fl2003h . We calculate an MSTO magnitude of M^gQy^, = 4.02 ±0.08 for NGC 



6752. 



We performed similar calc ulations using the NGC 0104 and NGC 6397 distance moduli 



listed in iGratton et al. I (120031 ) of (m - M)y = 13.50 ± 0.08 and 12.58 ± 0.08 , respectively. 
We calculate MSTO absolute magnitudes of M^gQy^, = 4.03±0.08 and M™g^^ = 3.77±0.08 
for NGC 0104 and NGC 6397, respectively. As the relative MSTOs are known, NGC 6752's 
MSTO absolute magnitude can be also derived independently from NGC 0104 and NGC 
6397 (3.72 ± 0.08 and 3.86 ± 0.08, respectively). The mean of the three NGC 6752 MSTO 
absolute magnitude determinations is MpQQQ^r = 3.87 ± 0.15 and we adopt this as the final 
MSTO zeropoint. 

We note that any error in the MSTO magnitude zeropoint would affect absolute age 
determination, but not relative age measurements. For this reason, the zeropoint uncertainty 
is not taken into consideration during relative age determination. 

Using the adopted zeropoint, we calculate MSTO absolute magnitudes for our entire 
sample, which are listed in column 4 of Table |H The listed uncertainties are determined in 
next section. 



5. Main sequence turn off uncertainties 



Several effects must be taken into account when deriving the uncertainty associated with 
the measured MSTO magnitudes {M^^q^). Among these, we address the cluster differential 
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reddening, the binary fraction and total number of cluster stars. The total number of stars 
affects the MSTO magnitude measurements in two ways. The first one is related to the total 
number of stars close to the MSTO that can be used to measure its magnitude: the smaller 
this number, the more uncertain the MSTO magnitude will be. On the other hand, a large 
number of cluster stars implies more severe crowding, and therefore less accurate magnitude 
measurement. To these uncertainties, we need to add the uncertainty from the MS— fitting 
and the mean ridge line determination procedure itself. The final uncertainty associated to 
the MSTO magnitude, ^^(M^gOgw), is: 

O (^F606w) ~ '^RED + ^BIN + ^NUM + '^CROW + ^MRL + '^MSF (2) 

where (T^ed the variance due to the differential reddening, a'^y^ due to the binary star 
population, cr^uM due to the number of cluster stars, ctq-^q^ due to the crowding, ct^rl is 
the MSTO's magitude variance from the mean ridge line determination, and ct^sf from the 
MS— fitting method. It is worth noting that field star contamination does not affect the 
MSTO determination. Our fields cover the central region in each cluster and the ratio of 
cluster stars to field stars is very high. Moreover, as noted above, the mean ridge line fitting 
method fits the dominant population, ignoring both field stars and blue stragglers. 

As it has been described in Section WA\ each MSTO's magnitude has been determined 
as the mean of the MSTO's magnitude obtained from ten different mean ridge line fittings. 
Here we assume that the variance associated with the mean ridge line determination, ct^rl, 
is the variance of the previous ten measurements. Typical values of ctmrl are in the range 
~0.01-0.04 mag. 

With the aim of performing an external cross-check for the ct^rl, the MSTO's magnitude 
variance from to the mean ridge line determination, we have considered four clusters having 
low reddening and a large number of stars (NGC 0104, NGC 5024, NGC 5272 and NGC7089). 
These clusters' CMDs have been randomly divided into 5 partial CMDs. A mean ridge line 
has been fitted to each partial CMD and its MSTO's magnitude has been measured. Finally, 
cr/-\/5 — 1 for these five MSTO's measurements has been computed for each cluster, where 
a is the standard deviation of the five measurements. 

Results are shown in Table [2l Column 1 shows the cluster's name, column 2 lists ctmrl 
and column 3 lists (t cross-check = cr/v5 — 1 derived in this external cross-check. It can be seen 
that the uncertainties obtained here are equivalent to the derived ctmrl- This test indicates 
that the adopted ctmrl is a good estimation. 

For cr^gp, we adopted the uncertainty calculated by the least squares main sequence 
fitting procedure. Typical values of ctmsf are in the interval ~0. 01-0. 06 mag. 
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To estimate the uncertainties induced by differential reddening, binary star population, 
and total number of cluster stars, as well as by differences between the measured MSTO 
and the actual GGC's MS TO, we generated over 250 synthetic CMDs using the lAC— star 



(lAparicio fc Gallart II2004I ) synthetic CMD program. Stellar populations are calculated on a 
star-by-star basis. L, Teg, and g of each star are computed by direct bilogarithmic interpola- 
tion in both the metallicity and age grids of a library of stellar evolution tracks, metallicity 
and age following continuous distributions. All the stellar evolution phases are covered accu- 
rate ly, and mass loss is compu ted during the RGB and the asymptotic giant branch phases 



see 



Aparicio fc Gallart I l2004j . for details). A set of different star formation rates, initial 
mass functions and metallicity laws are allowed, and binary star content can be computed. 
The result is a synthetic CMD with a smooth realistic stellar distribution. 

Synthetic CMDs were generated with an age of 13 Gyr, [M/H] = —1.2, differen- 
tial reddening E{F606W — F814W) values between 0.0 and 0.2, binary fractions (with 
0.5 < m2/mi < 1.0) between 0.0 and 0.7, and a number of stars brighter than Mpeoew = 6 
(approximately two mag. below the MSTO) between 1,000 and 50,000. The P04 stellar evo- 
lution library was used to generate the stellar libraries. For each combination of differential 
reddening, binary fraction and number of stars, 10 different CMDs were generated. Obser- 
vational errors were simulated in all the synthetic CMDs by applying a gaussian dispersion 
to the synthetic CMD, with a taken from the dispersion of the observed CMD as a function 
of the magnitude. 

Using the same mean ridge line fitting method used for the observational data, MSTO's 
magnitudes were measured for the synthetic CMD. The results are shown in Figure |H Grey 
points show the MSTO measurements from the 10 different CMDs generated for each exper- 
iment, black points show the mean, and the error bars correspond to the standard deviation 
of the 10 measurements. Solid line represents the input MSTO magnitude, that is, the 
magnitude of the MSTO obtained directly from the input isochrone. It can be seen that a 
systematic effect on the measured MSTO's magnitudes exists. Output MSTOs (measured) 
are systematically ~0.02 mag. fainter than input ones. This would have an effect on the 
measured absolute ages, but as this work is focused on relative ages, the variance of the 
difference between the input and output MSTO's magnitudes, cxjo, will be also considered. 

The left panel of Figure H] shows the Mpg^g-y^ measurements for the synthetic CMDs 
with no differential reddening, no binary stars and varying the number of stars brighter than 
-^F606W = 6 from 1,000 to 50,000. With the aim of determining cTnum; ^ gaussian has been 
generated for each black point, being the mean and the sigma of the gaussian equal to the 
value of the corresponding point and its sigma, respectively. Finally, all gaussians have been 
added together and the variance of this sum of gaussians has been adopted as ct^um- Since 
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most of the target clusters have more than 1,000 stars brighter than FQ06W = 6, we adopt 
c"num = 1-6 X 10"^ mag. (dotted hnes). Estimated this way, ct^um ^^^o includes the ctjq 
contribution. 

The central panel of Figure H] shows the effect of differential reddening upon the esti- 
mated MpgQg^ value. For a synthetic CMD, the reddening has been simulated by shifting 
every star along the reddening vector a quantity between zero and a maximum differential 
reddening value, E{F606W — F814:W)max- We generated a set of nine synthetic CMDs with 
50,000 stars brighter than F606W = 6 and different values of E{F606W - F8UW)max from 
0.0 to 0.2 magnitudes. It is worth noting that the black points' variance also includes the 
contribution from a^Q. In order to avoid duplicity, the contribution of a^g must be elimi- 
nated. As these synthetic diagrams have 50,000 stars brighter than FQ06W = 6, their ctj^q is 
equal to the variance of the point in the left panel corresponding to a number of stars brighter 
than F606W = 6 equal to 50,000, that is aio = 9.7 x 10"^ mag. In order to compute the 
effect of differential reddening upon the estimated Mpg^g^ value, a^Q has been quadratically 
subtracted to each point's variance in the central panel of Figure |H The resulting variance 
accounts for the differential reddening effect only. Then a gaussian has been generated for 
each black point, being the mean and the sigma of the gaussian equal to the value of the 
corresponding point and the obtained quadratically subtracted sigma, respectively. Finally, 
all gaussians have been added together and the variance of this sum of gaussians has been 
adopted as ct|^ed- we expect most of the target clusters to have differential reddening 
lower than 0.15 magnitudes, we adopt ctred = 2.0 x 10~^ mag. 

The right panel of Figure |4] shows the Mpg^g-y^ measurements of the synthetic CMDs 
with 50,000 stars brighter than Mpeoew = 6, no differential reddening, and varying fractions 
of binary stars. It can be seen that the fraction of binary stars affects Mpg^g^, making 
the measured MSTO fainter for clusters with a higher binary fractions. Following the same 
procedure as before, ctjq has been quadratically subtracted to each point's variance in the 
right panel of Figure HI Then a gaussian has been generated for each black point, being the 
mean and the sigma of the gaussian equal to the value of the corresponding point and the 
sigma obtained by subtracting (Tjq to the rms of the point, respectively. Finally, all gaussians 
have been added together and the variance of this sum of gaussians has been adopted as 
cr|jjv^. In a parallel work, we are measuring the fraction of binary stars in each cluster in our 



database (IMilone et al. Il2008bl ). For most clusters, it is smaller than 20%. Accordingly, we 
consider the fraction of binary stars interval [0.0, 0.2] which results in ctbin = 5.9 x 10^'^ 
mag. 



To estimate the effect of cluster crowding on the r neasu red MSTO magnitude, we used 
the artificial— star tests described in lAnderson et al. I (120081 ). In principle, crowding would 
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produce some systematic shift between output and input magnitudes of artificial stars. We 
are interested in evaluating if these shifts significantly change as a function of crowding. Since 
different clusters have different crowding levels, determining the artificial star magnitude 
shifts in several of them will provide the information we need. We will use the NGC 0104, 
NGC 1851, NGC 2808, NGC 5139, NGC 5286, NGC 6171 and E3 clusters. This sample 
contains both one of the most and one of the least crowded clusters, namely NGC 2808 and 
E3, respectively. 

Each cluster has associated a table containing input and output positions and magni- 
tudes of the artificial stars. For each cluster, we selected artificial stars with input F814W 
(/) magnitudes within ±0.25 mag of the MSTO. For most clusters, ~ 4, 000 artificial stars 
fulfill this condition. For each of them, we calculated the quantity 61 = lout — hni where 
lout and lin are the recovered and input magnitudes, respectively. To have a representative 
estimate of 51 and also of (j{51) for each cluster, we have divided the artificial stars in 10 
subsamples and calculated med(5/) for each of them. The average of medians, < med(5J) >, 
and the mean square root, (j[med(5/)], of the 10 med(5/) values have been computed for 
each cluster. Standard errors (T„ = cr[med((5/)]/ \J {n — 1), with n = 10, are obtained. These 
can be considered as an estimate of the internal errors of med(5/) of each cluster. Average 
values < med(5J) > run from 3.0 x 10~^ mag for E3, to —4.5 x 10~^ mag for NGC2808, while 
(Tn is always of the order of 1.0 x 10""^ mag, indicating that the med(5/) are significant. In 
summary, we can say that ~ 5 x 10~^ mag is the maximum shift introduced by crowding in 
our MSTO measurements. Finally, we can estimate the uncertainty produced by crowding 
that we were looking for, (Tcrow! as the root mean square of the < med(5/) > values of each 
cluster. This results in ctcrow = 1-5 x 10~^ mag. 

Table [3] lists a summary of the different contributions to MpgQg^ uncertainties. The 
MSTO's magnitude variances produced by the cluster number of stars, differential reddening, 
binary stars content and crowding, evaluated in this section using simulations, have been 
quadratically added and the resulting sigma (0.026 mag.) is also shown in Table [31 Finally, 
^F606W uncertainty has been evaluated for each cluster adding quadratically 0.026 mag. to 
the obtained ct^rl cr^gp (equation [2]), and the results are listed in Table H] (column 4). 
It is worth mentioning that the total uncertainties tend to increase both at low and high 
metaUicities. This is due to the error propagation in the cr^gp contribution, as we are using 
the intermediate metallicity cluster NGC 6752 for the MSTO absolute magnitude zeropoint. 
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6. Results 

In this section, we present the relative ages derived for the sample of 64 GGCs. We 
compare the derived ages to the theoretical isochrones of D07 using both the ZW and CG 
abundance scales. In the second part of this section, we test the MS— fitting procedure and 
check the self— consistency of the relative age determination. The relative age results are 
then compared with those derived using the stellar evolution libraries of P04, B94, and GOO. 
Finally, we present a comparison of our results with previous work. 



6.1. Relative ages 



Relative ages were calculated using a stellar evolution library (DOT) to calculate the 
theoretical M^^q^ values for different ages and [M/H]. Figure O shows the M^^q^ resulting 
from the DOT models as a function of [M/H], using both ZW and CG metallicity scales. Lines 
represent the DOT model MSTO magnitudes in steps of 1 Gyr (solid lines) and 0.5 Gyr (dotted 
lines) . The curves are interpolated using a spline curve so that we can easily estimate Mpg^g^ 
= f([M/H], age). We overplot the turnoff magnitudes calculated in section (open circles). 
The GGC metallicities listed i n Table [J have beeri transformed into global metallicities, 
[M/H], using the prescription of ISalaris et al. I (Il993l ): 



[M/H] = [Fe/H] + log{0.Q38f + 0.362) 



(3) 



where loaf = [a/Fe]. For the g-element enhancemerit, based on previous litera ture estimates 



(ICarney Ill996l : ISalaris fc Cassisi Ill996l : IVenn et al. II2004I : iKirby et al. II2008I ). we assumed: 
[a/Fe] = +0.3 ± 0.05 for clusters with [Fe/H] < -1.0, and [a/Fe] = +0.2 ± 0.05 for 
clusters with [Fe/H] > —1.0. The effect that this ± 0.05 dex uncertainty has on the final 
relative ages will be discussed at the end of this section. The values of [M/H] in both ZW 
and CG scales are listed in columns 2 and 3 in Table HJ respectively. We finally estimate the 
age of each cluster based on the interpolated curves. 

While this procedure provides absolute ages for all of our clusters, it must be noted 
that absolute ages depend on the theoretical model as well as the adopted MSTO magnitude 
zeropoint. A much more detailed analysis is required to derive reliable absolute ages, and 
this will be done in a forthcoming paper. This paper concentrates only on relative ages, 
which are much less dependent on the MSTO zero point. 



In order to derive relative ages, we divided the target clusters into three groups: low— metallicity 
([M/H]<-1.4), intermediate-metallicity (-1.4<[M/H]<-0.8), and high-metallicity ([M/H]>- 
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0.8). Using the D07 models, we compute the mean age of the low— metalhcity group as 
12.80 ± 0.17 Gyr. The rms scatter of the 13 low— metallicity GGCs is equal to 0.6 Gyr. The 
absolute age of each cluster was then divided by this mean age, obtaining what we will call 
normalized age from here on. Columns 5 and 6 of Table H] list the normalized ages from the 
D07 models using the ZW and CG metallicity scales, respectively. 

The adopted ±0.05 dex uncertainty for the a-element enhancement translates into an 
uncertainty of ~ 0.012 on the final relative age. With the aim of taking into consideration 
this effect on the relative ages, this quantity has been added quadratically to obtain the final 
relative age uncertainties listed in Table |H 

6.2. Testing self— consistency 

In this section, the self— consistency of the MS— fitting procedure and relative age de- 
termination is tested. In particular, we are interested in quantifying the validity of matching 
up clusters with different metallicities on the grounds of theoretical stellar evolution models. 

We created a set of 13 Gyr synthetic isochrones using the D07 models with metallicities 
similar to those of the selected reference clusters NGC 6304, NGC 6723, NGC 6981, NGC 
6681, NGC 6101 and NGC4590. That is, isochrones with [M/H] = -0.5, -0.8, -1.0, -1.1, -1.5 
and -1.8 have been considered. We applied the same MS— fitting procedure to the isochrones 
as for the observational data (section [4. 2p . The [M/H] = -1.1 isochrone, located in the center 
of the metallicity interval, was chosen for the MSTO zeropoint determination. 

Figure summarizes the test results. Upper panel shows the MS— fitting for the 13 Gyr 
synthetic isochrones with the same metallicities as the adopted reference clusters. Lower 
panel shows the obtained relative ages results. The MSTO of the [M/H] =-1.1 isochrone, 
that has been adopted as MSTO zeropoint, has been encircled. Error bars represent the 
uncertainties derived from the MS— fitting, (Tmsf- It can be seen that the main sequence 
fitting procedure is self consistent, and that the derived ages are similar to the input ages with 
a typical uncertainty of less than 2% over the entire metallicity range. We thus conclude that 
our metallicity grouping should not induce a significant bias in the relative age estimation. 

6.3. Comparison to other stellar evolution libraries 

Figure [7] shows the theoretical MSTO magnitudes derived using P04, B94, and GOO 
stellar evolution libraries. D07 model predictions have been also plotted for comparison. 
Lines represent Mpg°g^ for different ages as a function of global metallicity, in 1 Gyr (solid 
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lines) and 0.5 Gyr (dashed lines) age steps. Points represent our measured MSTO magnitudes 
as a function of metallicity in the CG scale. In the case of GOO models, metallicities lower 
than Z=0.0004 have been extrapolated. 

Using these theoretical grids, normalized ages have been derived following the procedure 
described in section 16. 1[ The corresponding average age of the low metallicity group, as 
derived from each isochrone set, has been used for normalization. 

Columns 7, 8 and 9 of Table H] list the normalized ages derived from the P04, B94 and 
GOO stellar evolution libraries, respectively, using the CG metallicity scale. Figure [8] shows 
the difference between the normalized ages from the D07 models, and those from the P04, 
B94 and GOO models, adopting the CG metalhcity scale. The error bars are the quadrature 
of the relative age's uncertainties form the two corresponding models. Interestingly, the 
relatives ages are the same, within the uncertainties, independent of the adopted theoretical 
library. 



6.4. Comparison to previous studies 



Figure M shows a comparison of our results to previous work. The upper and mid- 
dle panels show t he di fference between our results and the normalized ages published by 
De Angeli et al. I (120051 ) using HST snapshot -deA05S- and ground base d -deA05G- data. 



The l ower panel shows the difference between our results (MF08) and [Rosenberg et al. 



( Il999l )'s normahzed ages. T he uncertainties a re the qua drature of our listed unc ertainties 
and those published i n the iDe Angeli et al. I (120051 ) and [Rosenberg et al. I (119991 ) studies. 
The jPe Angeli et al. I (120051 ) ages are consistent with ours within the error bars and no ap- 
preciable age trend is observed. On the oth er hand, some trend is see n when comparing our 
results to those by iRosenberg et al. I (Il999l ): iRosenberg et al. I (Il999l ) find marginally older 



ages for cluster with —1.4 <[M/H]< —0.9 and younger ages for high metallicity clusters. 
This discrepancy is in part due to th e different methods app lied to measure the observed 
cluster parameters (see discussion in iDe Angeli et al. I l2005l . and in particular their Fig- 



ure 7), and in part to the different evolutionary models adopted for the age determination 
(IDe Angeli et al. II2005I . uses the P04 models, which we have shown to give consistent results 
with our adopted D07 models) . 
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6.5. The special case of NGC 0288 and NGC 0362 

In this study MSTO's magnitude is measured for each GGCs in our database, and then 
obtained magnitudes are transformed into ages using a set of theoretical stellar evolution 
models. During this transformation, canonical chemical abundances are assumed for all 
GGCs. We are aware of the fact that if a particular cluster has different chemical abundances, 
for example, an anomalous CNO content, then its age determination could be affected by 
this effect. It is worth mentioning that age uncertainties could be underestimated in this 
study because of the uncertainties coming from the chemical inhomogeneities which, to date, 
are impossible to quantify. 

In particular, our results indicate tha t NGC 0288 and NGC 0362 have the same age 



withi n ± 0.9 Gyr, while previous works (ISarajedini &: Demarque I Il990l : iBellazzini at al. 



20011) have found NGC 0288 to be 1-2 Gyr older than NGC 0362. The relative ages from 
these studies come from a comparison of all the evolutionary sequences (MS, SGB, HB, 
etc.), and not on a measurement of the absolute magnitude of the MSTO, as in the present 
paper. Therefore, we are planning a forthcoming paper in which we will, among other things, 
consider the case of these two clusters in more detail. In this context, it is worth mentioning 
that the present study does not take into consideration the distribution of stars along the 
HB, or finer features such as the RGB tip luminosity, the RGB bump, AGB bump and 
multiple populations. 

A similar argument could be used for the cas e of NGC 5272 and NGC 6205, for 



which we obtain identical age, while prev ious works (IVandenBerg. Bolte. &: Stetson I Il990 



Chaboyer. Demarque. &: Sarajedini Ill996l ) have found NGC 6205 to be 1-2 Gyr older than 



NGC 5272. 



7. Discussion 



Figure [To] shows the GGC normalized ages derived using the D07 stellar evolution library 
as a function of [M/H] in the CG metallicity scale (upper panel), and as a function of the 
galactocentric distance {tqc, lower panel). The low—, intermediate— and high— metallicity 
subgroups are separated by long dashed lines, and clusters belonging to the different metal- 
licity groups are plotted with open circles (low—), filled triangles (intermediate—) and filled 
circles (high— metallicity). The mean and dispersion of the low— metallicity clusters' age are 
marked with solid and dashed lines, respectively. For each of the three metallicity groups, 
mean age and rms is indicated. Figure [TT] shows the same data as Figure [TU| but for the ZW 
metallicity scale. 
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Overall, we notice an increase in age dispersion (which is model independent) with the 
metallicity. However, we can look at the results shown in the upper panels of Figures [10] 
and [TT] in another way. They show two branches in the age— metallicity relation of GGCs. 
On one hand, a 'you^^g' branch showing a clear age— metallicity relation, age decreasing at 
higher metallicities. On the other hand, a largely coeval 'old' branch. The 'old' branch shows 
an apparent relative age dispersion of ~0.05 and no age— metallicity relation. It is worth 
mentioning that while the age— metallicity relation may be model dependent (though it is 
somehow reassuring that the four most recent theoretical libraries provide consistent results 
on this respect) the age dispersion— metallicity relation is not model dependent. From now 
on, we will divide the clusters in two groups that we call 'young' and 'old'. GGCs with total 
metallicity higher than -1.4 and normalized age younger than 0.95, that is one sigma smaller 
than the mean of the low— metallicity clusters' age, are considered members of the 'young' 
group. The remaining GGCs are considered members of the 'old' group. The only exception 
is NGC 6441, a very high metallicity cluster which is much closer to the old branch than the 
young one, so it is included in the old group. Column 10 in Table H] lists whether the target 
cluster belongs to the old or young group. 

An age trend with galactocentric distance is seen in the lower panels of Figures [TD] 
and [m The fraction of young group's clusters increases significantly as the galactocentric 
distance increases. As a result, age dispersion increases at increasing galactocentric distance. 
Other authors ha ve commented on t he increased dispersion in ages at larger galactocentric 



distances, (e. g. [Richer et al. II1996I . which is one of the first to state the empirical problem 
clearly). Also in these panels, a clear metallicity trend with galactocentric distance can be 
observed for the old group of clusters, with an increasing metallicity at lower galactocentric 
distances. This trend is explicitly shown in Figure [T2l which plots GGCs [M/H] in the CG 
metallicity scale as a function of the galactocentric distance. Open and filled circles represent 
clusters in the young and old groups, respectively. 

In the following, the age dispersion of the old clusters group, the age— metallicity relation 
of the young clusters group and the age trend with galactocentric distance are discussed in 
detail. 



7.1. The age dispersion of the old clusters group 

Figures [10] and [11] indicate that the mean normalized age of the clusters in the old group 
is the same for all metallicities. The relative age dispersion of this group is ~0.05. If an 
absolute age of 12.8 Gyr (from the D07 models) is assumed as the mean age of these clusters, 
age dispersion would be ~0.6 Gyr. We note that if we assume that all sources of uncertainty 
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have been taken into account in the error bars, then at least part of this dispersion must be 
real. The average normalized age uncertainty among the clusters in the old group is 0.04, 
which is approximately 0.5 Gyr. Subtracting this quadratically from the observed dispersion, 
we find an intrinsic dispersion of 0.03 in relative age for the clusters in the old group. Using 
12.8 Gyr for the age of this group, the intrinsic dispersion is ~0.4 Gyr. 

These results are consistent with a scenario in which the old group of clusters formed 
within a fast assembling process of the halo, lasting ~0.8 Gyr and which should also ac- 
count for an increasing metallicity towards the center of the Galaxy. This implies that the 
chemical enrichment of the protogalactic cloud was faster than the assembling timescale. 
Regarding the natu r e of s uch assembling process, it is interesting to note the following. 
Wilkinson &: Evans I (119991 ) estimate the mass and scalelength of the Milky Way dark mat- 
ter halo to be ~ 1.9 x 10^^ Mq and 170 kpc, respectively. The free-fall time of a homogeneous 
sphere of those mass and radius is ~ 0.84 Gyr, similar to the age range we find here for the 
old clusters. In other words, the age dispersion of old globular clusters is not in contradiction 
wi th the formation from the colapse of a si ngle protosystem, resembling the model proposed 
by lEggen. Lynden-Bell. fc Sandage I (1l962l ). On another side, the standard ACDM scenario 
would predict a general star formation, including star clust ers, in protogalaxy building blocks 
of mass ~ 10^ Mq before reionization (IMoore et al. II2006I ). These blocks merge together af- 
terwards to form larger galaxies, and would, in principle, be compatible with the old, coeval 
cluster population that we observe. More difficult seems to account for the metallicity gra- 
dient (Fig. [T^ under this scenario. A more detailed analysis is necessary, but it is important 
to note that any successful galaxy formation scenario must account for the existence of a 
large nubmer of old, coeval globular clusters in the Milky Way. 



7.2. The Age— metallicity relation of the young clusters group 



Clusters in the young group, however, show a clear age— metallicity relationship, with 
younger clusters being more metal rich than older ones. The nature of this "young branch" 
can perhaps be better understood by looking at the GGCs associated with known or putative 
accreted dwarf galaxies: Sagittarius, Canis Major and Monoceros. 



Ibata. Gilmore. fc Irwin I (119951 ) discovered t he accreted Sagittariu s dw arf galaxy and its 



GC sy stem (Terzan 7, Terzan 8, Arp 2, and M54). iDinescu et al. I (120001 ) and lBellazzini et al. 
(120031 ) have since then argued, respectively, that Pal 12 and NGC 4147 are part of the ex- 
tended tidal stream. 



The Monoceros Ring was discovered by 



Newberg et al, 



(I2OO2I ). and it has been pro- 
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posed to be the t i dal st ream of the accreted Monoceros galaxy. ICrane et al. I (120031 ) and 
Frinchaboy et al. I (120041 ) identified five possible GC candidates of this accrete d gala xy: NGC 



2298, NGC 2808, NGC 5286, Pal 1, and BH 176. In addition. Martin et al. I (120041 ) identify 



four GGCs - NGC 1851, NGC 1904, NGC 2298 and NGC 2808 - possibly associated with 
a stellar overdensity in Canis Major, which they suspect to be a dwarf galaxy. However, 
the nature of these two structures is cont roversial. The existenc e of th e Ca nis Major dwarf 



galaxy has been seriously questioned by iMomany et al. I (120041 . |2006| ). In iMomany et al. 



( l2006l ). the Canis Major stellar over— density has been completely accounted for as the effect 
of the Milky Way disk warp (an ubiquitous property of all massive galaxies). In addition, 
Momany et al. I (120061 ) have shown that also the Monoceros Ring seems to refiect the signa- 
ture of another Galactic disk property; the fiaring (the increase in scale-height as a function 
of Galactocentric distance) of the outer disk. 

Figures [T31A andfTSB show the age— metallicity relation, using the CG metallicity scal^. 
Sagittarius, Monoceros and Canis Major's clusters have been marked. It is interesting that, 
for metallicities higher than [M/H] = -1.4, all of Sagittarius, Monoceros and Canis Major's 
accreted clusters fall in the young group, with the exception of NGC 5286, and they follow 
the same age— metallicity relation as the rest of the young group's clusters. This result 
suggests a different origin for the old and young groups of clusters. It is worth mentioning 
that there is a number of clusters, other than the Sagittarius, Monoceros and Canis Major's 
clusters among the young sample. 

Another interesting question is related to the recent discovery that GGCs are not sim- 
ple single stellar populations. The high mass and multiple stellar populations of uCen 
has led to speculation that it is the remnant nucleus of an accreted Milky Way satellite 
galaxy (Lee et al. [l999 : Villanova et al. 2007 ) . But c uCen is not a unique case. NGC 2808 



(IPiotto et al. 



20071). NGC 1851 (IMilone et al. Il2008ah . NGC 6388 (IPiotto Il2008f ) and NGC 



tions. 



Freeman 



6441 (ICaloi fc D'Antona 1 120071 ) are four massive clusters hosting multiple stellar popula- 



Jl993h suggested that GGCs form as the nuclei of dwarf galaxies in the 
early universe, and are accreted as their host galaxies merge onto larger structures. In this 
context, these peculiar GGCs would be good candidates to be the remnants of accreted 
satellites. Multiple population clusters have been also marked in Figure [T3JA.. 

We note that multiple stellar population clusters have broadened MSTOs and at least 
a fraction of the stars in these clusters have rather anomalous chemical compositions. For 
this reason, we should be hesitant in adopting their relative ages obtained in this study. A 



^Because the results are equivalent to those using the ZW metallicity scale, we only consider the CG scale 
after this point. 
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muc h more detailed ana lysis is needed for these peculiar clusters. iVillanova et al. I (120071 ) 
and ICassisi et al. I (120081 ) present particularly illuminating discussions of the difficulties in 



estimating ages for the multiple populations in cuCen and NGC 1851, respectively. 

Figure [T3B shows the same Figure [T31A. but here GGCs with tqc > 10 kpc have are 
represented with filled circles, while those with tqc < 10 kpc are plotted with open circles. 
Interestingly, all intermediate— and high— metallicity clusters with tqc > 10 kpc fall on the 
young group of clusters. 

In order to characterize the young group's age— metallicity relation, a least squares fit 
has been performed, and results are also shown in Figure [T3B. The age— metallicity relation 
observed in the young group of clusters can be described by the equation: 



AgeNORM = -0.38[M/H]+0.45 (4) 

This relation is plotted with a solid line in Figure [T3lB. The dispersion of the clusters' 
relative ages with respect to the previous equation is also 0.05, and it is also represented in 
Figure [TSB (dashed lines). Following a similar argument as for the old group, if an absolute 
age of 12.8 Gyr is assumed as the mean of the clusters in the old group, the young group's 
age dispersion with respect to the previous age— metallicity relation would be ~0.6 Gyr. The 
average normalized age uncertainty among the clusters in the young group is 0.04, which is 
approximately 0.5 Gyr. Subtracting this quadratically from the observed dispersion, we find 
an intrinsic dispersion with respect to its age— metallicity relation of 0.03 in relative age, or 
~0.4 Gyr. 

In summary, our results are consistent with the clusters of the young group having been 
formed in a different process spanning an interval of time as long as 0.45 in relative age, or 
~6 Gyr, and resulting in a group of GCs with a clear age— metallicity relation. It is very 
tempting to argue that the origin of this second group of clusters is related to their formation 
within Milky Way satellite galaxies that were later accreted. However, the reason why the 
age— metallicity relation of the GGCs associated with the Sagittarius dwarf galaxy, as well 
as Monoceros and Canis Major, is similar to that of the rest of 'accreted' clusters should 
be investigated. Perhaps all of the young group clusters share the same origin, which would 
require an association between Sagittarius, Monoceros and Canis Major, or alternatively, 
all dwarf galaxy systems may share a common (or very similar) star formation and metal 
enrichment history, and hence a common age— metallicity relation, which seems unlikely. 
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7.3. Age trend with galactocentric distance 

With the aim of analyzing in more detail the observed age trend with galactocentric 
distance, GGCs in the old and young groups are considered separately. Figure [T3D shows 
normalized ages versus galactocentric distance (same as Figure [TOl lower panel), but now 
clusters in the old group are represented with open circles and clusters in the young group 
with filled circles. Besides, Sagittarius, Monoceros and Canis Major's clusters are shown as 
encircled points. 

It can be seen that if the old group is considered, an age trend with galactocentric dis- 
tance is not observed. The age's variance remains constant with the galactocentric distance 
for this group. 

On the other hand if the young group is taken into account, it can be seen that the 
age's variance increases with galactocentric distance. In order to determine if this variance 
increase depends either on the galactocentric distance or the metallicity, a three dimensional 
principal— component (PC) analysis has been done considering the clusters in the young 
group only. The obtained eigenvalues are 2.03, 0.83 and 0.14. The first PC, of eigenvalue 
2.03, accounts for 68% of the total variance. According to the commonly used average or 
eigenvalue-one criterion, only this PC should be retained. Figure [T3] shows the relations 
between the original data and PCs projected on the two first PCs plane. On the one hand 
it can be seen that the first PC is strongly correlated with [M/H] and age, but weakly with 
tgci while the second PC carries on most of the roc variance. On another hand it is clear 
that young cluster age dispersion can be explained by an age— metallicity relation and that 
age is not related to tgc- 

In fact, the effect of metallicity on age of clusters in the young group can be eliminated 
by subtracting the fitted age— metallicity relation (equation H]) to the cluster's ages. By 
doing this, age— metallicity 'corrected' ages are obtained. Figure [T3D shows the same as 
Figure \WC but here age— metallicity 'corrected' clusters are represented with solid circles. 
It can be clearly seen that most of the age dispersion observed in Figure WSC disappears, 
and that age is not related to tgc 

In summary, it can be concluded that a significant part of the age dispersion present in 
Figure [T3P for tgc ^ 15 kpc and most of that for tgc > 15 can be explained by a strong 
age— metallicity relation for the clusters of the young group. 

It is worth mentioning that our database is biased in distance. No outer— halo clusters, 
with galactocentric distances larger than ~20 kpc, have been considered for this study. So 
the analysis of the galactic halo's age structure is limited by the galactocentric distance range 
of the database. Given that outer— halo clusters (Pal 3, Pal 4, Pal 14, Eridanus, AM-1) have 
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younger ages (ISarajedini 1119971 : IStetson et al. Ill999l : iDotter et al. II2008I ). it is possible that 



our age range might be a lower limit on the actual range. 



8. Conclusions 

Normalized ages have been derived for a sample of 64 GGCs using the stellar evolution 
models of DOT and both ZW and CG metallicity scales. We have also performed an analysis 
using the stellar models of P04, B94 and GOO. The result is the most extensive and precise 
database of normalized ages so far produced . Our results are: 

• We find that we are able to measure relative ages to a formal precision of 2-7%. The 
relative ages are independent of the choice of theoretical model. Four independent sets 
of stellar evolution libraries (D07, P04, B94, GOO) produce essentially identical results. 

• We find that the GGCs fall into two well-defined groups. The first one represents 
a population of old clusters that have the same age, a dispersion of 5% in relative 
age and no age— metallicity relation. If we assume an absolute age of 12.8 Gyr, the 
absolute age dispersion of the old clusters is ~0.6 Gyr. Accounting for the measurement 
uncertainties, we obtain an intrinsic age dispersion of ~0.4 Gyr. The second group of 
clusters shows a clear age— metallicity relation, with young clusters being more metal 
rich than older ones. Also in this case, the intrinsic age dispersion with respect to its 
age— metallicity relation is ~0.4 Gyr. Besides, there is no age dispersion— metallicity 
relation if we look at the two samples (old and young) separately. 

• These results are consistent with a scenario in which the formation of the Milky Way 
CCG system took place in two phases. The first one produced clusters with ages within 
a range of ~0.8 Gyr. This age range is compatible with the timescale for the collapse of 
a protogalaxy of the same mass and scalelength as that of the Milky Way dark matter 
halo. In other words, the age dispersion of the old group of globular clusters is not in 
contradiction with the formation from the colapse of a s i ngle p rotosystem, resembling 



the model proposed by lEggen. Lynden-Bell. fc Sandage I (119621 ). The standard ACDM 



scenario would also predict a significant star formation in protogalactic building blocks 
before reionization, which, after merging, would produce an old, coeval population of 
globular clusters in a large galaxy like the Milky Way. However, to account for the 
galactocentric metallicity gradient of this group seems difficult in the context of this 
mechanism. The second phase spanned a time interval as long as ~6 Gyr and resulted 
in a group of GGCs with a clear age— metallicity relation. It is very tempting to argue 
that the origin of this second group of clusters is related to their formation within 
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Milky Way satellite galaxies that were later accreted. However, this would not account 
for the fact that the clusters of this group share the same age— metallicity relation. 

• We find that the increasing dispersion in age with galactocentric distance is explained 
by the age— metallicity relation of GGCs. 

• These results are independent of the assumed metallicity scale. 

• It is worth mentioning that our database is biased by distance. No outer halo clusters, 
with galactocentric distances larger than ~20 kpc, have been considered. The analysis 
of the galactic halo's age structure is therefore limited by the galactocentric distance 
range of the database. 
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(F606W-F814W) 



Fig. 1. — Mean ridge line fitting procedure using rototec? histograms. This shows the analysis 
of NGC 1851. In short, the mean ridge line is determined in two steps. First, color histograms 
are computed for stars grouped in 0.4 magnitude (F8141^) wide bins (left panel). A prehm- 
inary mean ridge line is then derived with {F606W — F814W) colors set at each histogram's 
maximum, and the F81AW magnitude of the bin center. Second, we use rotated histograms, 
perpendicular to the preliminary mean ridge line, to derive the final one (right panel). See 
text for details. 
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Fig. 2. — Examples of mean ridge line fitting in four selected clusters covering the metallicity 
range. The MSTO position is marked with an open circle. 
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(F606W-F814W) 

Fig. 3. — Examples of MS— fitting for the —1.3 < [Fe/ H]cg < -1-1 metallicity group. The 
reference cluster is NGC 6981 (dashed hue). Each cluster MRL (solid hne) has been fitted to 
the reference cluster in the magnitude intervals [(MpgQg^^ — 2.5) < Mpeoew < (^Feoew ~ l-^)] 
and [(Mje^ew + 1-5) < Mpeoew < (Mpeoew + 3-0)] (shaded regions). 
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Fig. 4. — The right panel shows MjgQg^ measurements over synthetic CMDs with 50,000 
stars brighter than Mpgoew — 6, no differential reddening, and varying the fraction of binary 
stars. Grey points show the MSTO measurements for the 10 different CMDs generated 
in each case, and black points show the average and standard deviation of each of the 10 
measurement sets. Solid line represents the MSTO's magnitude measured directly on the 
isochrone. The Central panel shows the same as left panel, but now the fraction of binary 
stars is 0.0 and the maximum differential reddening is varying between 0.0 and 0.2. The left 
panel shows the same, but now the fraction of binary stars is 0.0, the differential reddening 
is 0.0 and the number of stars brighter then Mf606W = 6 is varying between 1,000 to 50,000. 
Dashed hues represent the adopted ±a in each case. See text for details. 
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Fig. 5. — Theoretical TO absolute magnitudes from the D07 stellar evolution library. Lines 
represent Mjgoew different ages as a function of global metallicity, in 1 Gyr (solid lines) 
and 0.5 Gyr (dashed hues) age steps. The lower isochrone corresponds to an age of 15 
Gyr. Points represent the measured MSTO magnitudes. The upper panel shows the Mjggg^ 
magnitude as a function of the metallicity in the ZW scale; the lower panel gives Mp^g^ as 
a function of the metaUicity in the CG scale. 
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Fig. 6. — Testing the used main sequence fitting procedure and relative ages determination 
using 13 Gyr synthetic isochrones. Upper panel shows the MS— fitting for 13 Gyr synthetic 
isochrones with the same metallicities as the adopted reference clusters ([M/H] = -0.5, -0.8, - 
1.0, -1.1, -1.5 and -1.8). Lower panel shows the obtained relative ages. The MSTO magnitude 
of the [M/H]=-l.l isochrone (encircled point) has been considered as MSTO zeropoint. 
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Fig. 7.— MSTO magnitudes predicted by the D07, P04, B94, and GOO stellar evolution 

libraries. Lines represent M'^^q^/^, for different ages as a function of global metallicity, in 1 
Gyr (solid lines) and 0.5 Gyr (dashed lines) age steps. The lower isochrone corresponds to 
an age of 15 Gyr and the upper one to 6 Gyr. In the case of GOO models, metallicities lower 
than Z=0.0004 have been extrapolated. Points represent observational data. 
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Fig. 8. — Comparison of normalized ages derived from different sets of models in the CG 
metallicity scale. The difference between normalized ages derived using P04, B94 and GOO 
models, respectively, and the D07 models are shown. It can be seen that relative ages derived 
using the main sequence fitting method are independent from the adopted model. 



- 39 - 



0.4 
0.2 


-0.2 
-0.4 



-i 1 1 r- 



-1 1 1 [- 



~\ I I r~ 



Mean = -0.04 ±0.07 



_l I I L_ 



_l I I L_ 



_l I I L_ 



n 1 1 T" 



1-^ I 



Mean = -0.02 ±0.08 



-1 1 1 r- 



-1 1 1 r 



0.4 
0.2 


-0.2 
-0.4 




_l I I L_ 



_l I I L_ 



_l I I L_ 



0.4 
0.2 


-0.2 
-0.4 



n 1 1 T" 



? I 



Mean = 0.00 ±0.11 



_l I I L_ 



"T 1 1 r- 



LI 



_l I I L_ 



n 1 1 T" 



■1.5 



[M/H] 



CG 



_l I I L_ 



-0.5 



Fig. 9. — Comparison with pre vious works using the C G metaUicity scale. The differences 
between the normahzed ages of iDe Angeh et al. jusing HST snapshot -deAOSS- and 

ground based -deAOSG- data) and Rosenberg et al. (jl999l ) (R99) and the results derived in 
this paper by using D07 models (MF08) are indicated. 
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Fig. 10. — Globular cluster normalized ages as a function of [M/H] in the CG mctallicity scale 
(upper pannel), and as a function of galactocentric distance (rcc kpc, lower pannel). These 
results have been derived using the D07 stellar evolution library. Open circles, filled triangles 
and filled circles represent GCs within the low—, intermediate— and high— metallicity groups, 
respectively. For each of the three metallicity groups, mean age and rms is indicated. See 
text for details. 
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Fig. 11. — Same as Figure fTOl but for the ZW metallicity scale. 
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Fig. 12. — Globular cluster's [M/H] in the CG metallicity scale as a function of the galac- 
tocentric distance. Open and filled circles represent clusters in the young and old groups, 
respectively. 
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Fig. 13. — Age— metallicity relation using the CG metallicity scale (A and B) and normalized 
ages versus galactoccntric distance (C and D). In A, clusters associated with Sagittarius, 
Monoceros and Canis Major have been marked, together with multiple stellar population 
clusters. In B, GGCs with tqc < 10 kpc are represented with open circles, while those 
with tqc > 10 kpc are plotted with filled circles. A least squares fit to the young group's 
age— metallicity relation is also plotted. In C, Sagittarius, Monoceros and Canis Major' 
clusters are also marked. Finally, D shows the same as C, but here ages of clusters in the 
young group have been age— metallicity 'corrected' using the age— metallicity relation shown 
in B. See text for details. 
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Fig. 14. — Young group principal— component (PC) analysis: Relations between original 
data and principal components on the first two principal components plane. It can be seen 
that the first PC is strongly correlated with [M/H] and age, but weakly with tgc-, while the 
second PC carries on most of the tqc variance. See text for details. 
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Table 1. Globular cluster parameters 



Ref. cluster Metallicity group 


Name 


[Fe/H]zw 


[Fe/H]cG 


TGC 


NGC 6304 -0.8<[Fe/HJcG <-0.3 


NGC 0104 


-0.71 


-0.78 


7.4 




NGC 5927 


-0.32 


-0.64 


4.5 




NGC 6304 


-0.38 


-0.66 


2.2 




NGC 6352 


-0.50 


-0.70 


3.3 




NGC 6366 


-0.58 


-0.73 


5.0 




NGC 6388 


-0.74 


-0.77 


3.2 




NGC 6441 


-0.59 


-0.60 


3.9 




NGC 6496 


-0.50 


-0.70 


4.3 




JNGC 6624 


-0.50 


-0.70 


1.2 




IN CjrL/ DUO ( 


-0.72 


-0.78 


1.9 




iNCjrL/ DOOO 


-0.58 


-0.73 


b. / 




LYJNLrAT 


-0.62 


-0.64 


4.2 




PAL 1 


-0.60 


-0.70 


17.0 






-0.05 


-0.56 


16.0 


NGC 6723 -l.l<[Fe/H]cG <-0.8 


NGC 0362 


-1.33 


-1.09 


9.4 




NGC 1261 


-1.32 


-1.08 


18.2 




NGC 1851 


-1.23 


-1.03 


16.7 




NGC 6121 


-1.27 


-1.05 


5.9 




NGC 6171 


-1.09 


-0.95 


3.3 




INVjrU bOOZ 


-1.18 


-0.99 


5.1 




NGC DD02 


-0.99 


-0.97 


2.8 




NGC 6717 


-1.33 


-1.09 


2.4 




NGC 6723 


-1.12 


-0.96 


2.6 




E3 


-0.80 


-0.83 


7.6 




PAL 12 


-0.82 


-0.83 


15.0 


NGC 6981 -1.3<[Fe/HJcG <-l-l 


NGC 0288 


-1.40 


-1.14 


12.0 




NGC 2808 


-1.36 


-1.11 


11.1 




NGC 3201 


-1.53 


-1.24 


8.9 




NGC 5904 


-1.38 


-1.12 


6.2 




NGC 6218 


-1.40 


-1.14 


4.5 




NGC 6254 


-1.55 


-1.25 


4.6 




NGC 6584 


-1.51 


-1.30 


7.0 




NGC 6715 


-1.54 


-1.25 


19.2 




NGC 6752 


-1.54 


-1.24 


5.2 




NGC 6981 


-1.50 


-1.21 


12.9 


NGC 6681 -1.5<[Fe/H]cG <-l-3 


NGC 4147 


-1.77 


-1.50 


21.3 




NGC 5139 


-1.59 


-1.35 


6.4 




NGC 5272 


-1.66 


-1.34 


12.2 




NGC 5286 


-1.70 


-1.41 


8.4 




NGC 5986 


-1.65 


-1.35 


4.8 




NGC 6093 


-1.75 


-1.47 


3.8 




NGC 6205 


-1.63 


-1.33 


8.7 




NGC 6535 


-1.78 


-1.51 


3.9 




NGC 6656 


-1.75 


-1.49 


4.9 




NGC 6681 


-1.64 


-1.35 


2.1 




NGC 6934 


-1.54 


-1.32 


12.8 
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Table 1 — Continued 



Ref. cluster Metallicity group 


Name 


[Fe/H]zw 


[Fe/H]cG 


roc 




NGC 7089 


-1.61 


-1.31 


10.4 




ARP 2 


-1.74 


-1.45 


21.4 


NGC 6101 -1.8<[Fe/H]cG <-l-5 


NGC 2298 


-1.91 


-1.71 


15.7 




NGC 4833 


-1.92 


-1.71 


7.0 




NGC 6101 


-1.95 


-1.76 


11.1 




NGC 6144 


-1.81 


-1.56 


2.6 




NGC 6397 


-1.94 


-1.76 


6.0 




NGC 6541 


-1.79 


-1.53 


2.2 




NGC 6809 


-1.80 


-1.54 


.3.9 




TERZAN8 


-1.97 


-1.80 


19.1 


NGC 4590 -2.3<[Fe/H]cG <-l-8 


NGC 4590 


-2.11 


-2.00 


10.1 




NGC 5024 


-2.04 


-1.86 


18.3 




NGC 5053 


-2.10 


-1.98 


16.9 




NGC 5466 


-2.22 


-2.20 


16.2 




NGC 6341 


-2.24 


-2.16 


9.6 




NGC 6779 


-2.20 


-2.00 


9.7 




NGC 7078 


-2.13 


-2.02 


10.4 




NGC 7099 


-2.05 


-1.92 


7.1 
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Table 2. External cross-check for the Umrt determination 



Name 


fMRL 


cross — check 


NGC 0104 


0.011 


0.011 


NGC 5024 


0.024 


0.013 


NGC 5272 


0.011 


0.012 


NGC 7089 


0.018 


0.006 
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Table 3. Summary of Mj^g^ uncertainty contributions 



Contribution 


mag. 




0"MRL 


~0.01-0.04 




cmsf 


~0.01-0.06 




cnum 


0.016 




cred 


0.020 


Quadratic sum 


cbin 


0.0059 


0.026 mag. 


fCROW 


0.0015 





Table 4. Globular cluster relative ages 



Name 


[M/H]zw 


[M/H]cG 


■""Feosw 


D07zw 


D07cG 


P04cG 


B94cG 


GOOcG 


Group 



NGC 


0104 


-0, 


,57 


-0, 


,64 


4.13 


± 


0.07 


0, 


,96 


± 


0.07 


1, 


,02 


± 


0.07 


1, 


,07 


± 


0.08 


1.02 


± 


0.07 


1, 


,05 


± 


0.09 


old 


NGC 


0288 


-1, 


,18 


-0, 


,92 


3.82 


± 


0.04 


0, 


,88 


± 


0.04 


0, 


,83 


± 


0.03 


0, 


,83 


± 


0.03 


0.82 


± 


0.03 


0, 


,83 




0.03 


young 


NGC 


0362 


-1, 


,11 


-0, 


,87 


3.81 


± 


0.05 


0, 


,84 


± 


0.04 


0, 


,81 


± 


0.04 


0, 


,80 


± 


0.04 


0.81 




0.04 


0, 


,81 


± 


0.05 


young 


NGC 


1261 


-1, 


,10 


-0 


.86 


3.80 


± 


0.05 


0, 


.83 


± 


0.04 


0, 


.80 


± 


0.04 


0, 


,79 


± 


0.04 


0.80 


± 


0.04 





.79 


± 


0.05 


young 


NGC 


1851 


-1, 


,01 


-0 


.81 


3.79 


± 


0.05 





.80 


± 


0.04 





.78 


± 


0.04 


0, 


,77 


± 


0.04 


0.78 


± 


0.04 





.75 


± 


0.04 


young 


NGC 


2298 


-1, 


,69 


-1, 


.49 


3.80 


± 


0.05 


1, 


.01 


± 


0.05 


0, 


.99 


± 


0.05 


0, 


,99 


± 


0.05 


0.99 


± 


0.05 


0, 


.99 




0.05 


old 


NGC 


2808 


-1, 


,14 


-0 


.89 


3.86 


± 


0.03 





.91 


± 


0.03 





.85 


± 


0.03 


0, 


,85 


± 


0.03 


0.85 


± 


0.03 





.85 


± 


0.02 


young 


NGC 


3201 


-1, 


,31 


-1, 


.02 


3.75 


± 


0.04 


0, 


,86 


± 


0.03 


0, 


,80 


± 


0.03 


0, 


,79 




0.03 


0.80 


± 


0.02 


0, 


.81 


± 


0.03 


young 


NGC 


4147 


-1, 


,55 


-1, 


.28 


3.74 


± 


0.04 


0, 


,93 


± 


0.04 


0, 


,89 


± 


0.04 


0, 


,88 


± 


0.04 


0.89 




0.04 


0, 


,89 


± 


0.03 


young 


NGC 


4590 


-1, 


,89 


-1, 


.78 


3.61 


± 


0.05 


0, 


,89 


± 


0.03 


0, 


,90 


± 


0.04 


0, 


,89 


± 


0.04 


0.90 


± 


0.04 


0, 


,91 




0.04 


old 


NGC 


4833 


-1, 


,70 


-1, 


.49 


3.80 


± 


0.05 


1, 


.01 


± 


0.05 


0, 


.98 


± 


0.05 


0, 


,98 


± 


0.04 


0.98 


± 


0.05 


0, 


.98 


± 


0.05 


old 


NGC 


5024 


-1, 


,82 


-1, 


.64 


3.76 


± 


0.05 


1, 


.01 


± 


0.05 


0, 


.99 


± 


0.05 


0, 


,99 


± 


0.04 


1.00 


± 


0.05 





.99 


± 


0.05 


old 


NGC 


5053 


-1, 


,88 


-1, 


.76 


3.69 


± 


0.05 





.95 


± 


0.04 





.96 


± 


0.04 


0, 


,97 


± 


0.04 


0.96 


± 


0.05 


0, 


.96 


± 


0.04 


old 


NGC 


5139 


-1, 


,37 


-1, 


.13 


3.82 


± 


0.05 


0, 


.94 


± 


0.05 


0, 


.90 


± 


0.05 


0, 


,89 


± 


0.05 


0.90 


± 


0.05 


0, 


.89 


± 


0.05 


young 


NGC 


5272 


-1, 


,44 


-1, 


.12 


3.83 


± 


0.04 


0, 


,96 


± 


0.04 


0, 


,89 


± 


0.04 


0, 


,88 


± 


0.04 


0.90 


± 


0.03 


0, 


.89 


± 


0.04 


young 


NGC 


5286 


-1, 


,48 


-1, 


.19 


3.89 


± 


0.05 


1, 


.04 


± 


0.04 


0, 


.98 


± 


0.04 


0, 


,98 


± 


0.05 


0.97 


± 


0.05 


0, 


.98 


± 


0.05 


old 


NGC 


5466 


-2, 


,00 


-1, 


.98 


3.74 


± 


0.05 


1, 


.03 


± 


0.05 


1, 


.06 


± 


0.05 


1, 


,06 


± 


0.06 


1.08 


± 


0.03 


1, 


.07 




0.05 


old 


NGC 


5904 


-1, 


,16 


-0 


.90 


3.83 


± 


0.03 





.88 


± 


0.03 





.83 


± 


0.03 


0, 


,83 


± 


0.03 


0.83 


± 


0.03 





.83 


± 


0.02 


young 


NGC 


5927 


-0, 


,18 


-0 


.50 


4.15 


± 


0.07 





.88 


± 


0.06 





.99 


± 


0.07 


1, 


,03 


± 


0.08 


0.98 


± 


0.07 


1, 


.01 


± 


0.11 


old 


NGC 


5986 


-1, 


,43 


-1, 


.13 


3.88 


± 


0.05 


1, 


,01 


± 


0.04 


0, 


,95 


± 


0.04 


0, 


,94 


± 


0.05 


0.93 


± 


0.04 


0, 


,95 


± 


0.05 


old 


NGC 


6093 


-1, 


,53 


-1, 


.25 


3.87 


± 


0.05 


1, 


,03 


± 


0.04 


0, 


,98 


± 


0.04 


0, 


,98 


± 


0.05 


0.97 


± 


0.05 


0, 


,98 




0.05 


old 


NGC 


6101 


-1, 


,73 


-1, 


.54 


3.78 


± 


0.04 


1, 


,00 


± 


0.04 


0, 


,98 


± 


0.04 


0, 


,98 


± 


0.04 


0.98 


± 


0.04 


0, 


,98 


± 


0.04 


old 


NGC 


6121 


-1, 


,05 


-0 


.83 


4.01 


± 


0.05 


1, 


.02 


± 


0.05 


0, 


.98 


± 


0.05 


0, 


,99 


± 


0.07 


0.97 


± 


0.05 


0, 


.98 


± 


0.05 


old 


NGC 


6144 


-1, 


,71 


-1, 


.52 


3.88 


± 


0.05 


1, 


.10 


± 


0.05 


1, 


.08 


± 


0.05 


1, 


,07 


± 


0.06 


1.05 


± 


0.05 


1, 


.07 


± 


0.05 


old 


NGC 


6171 


-0, 


,87 


-0, 


.81 


4.14 


± 


0.05 


1, 


.08 


± 


0.05 


1, 


.09 


± 


0.06 


1, 


,11 


± 


0.06 


1.07 


± 


0.06 


1, 


.13 




0.08 


old 


NGC 


6205 


-1, 


,41 


-1, 


.11 


3.84 


± 


0.04 





.97 


± 


0.04 





.91 


± 


0.04 


0, 


,90 


± 


0.04 


0.90 


± 


0.02 





.90 


± 


0.04 


young 


NGC 


6218 


-1, 


,18 


-0, 


.92 


3.99 


± 


0.04 


1, 


,05 


± 


0.04 


0, 


,99 


± 


0.03 


1, 


,00 


± 


0.05 


0.99 


± 


0.04 


0, 


,99 




0.03 


old 


NGC 


6254 


-1, 


,33 


-1, 


.03 


3.86 


± 


0.04 


0, 


,96 


± 


0.04 


0, 


,89 


± 


0.04 


0, 


,88 


± 


0.04 


0.90 


± 


0.03 


0, 


,88 


± 


0.04 


young 


NGC 


6304 


-0, 


,24 


-0, 


.52 


4.19 


± 


0.07 


0, 


,93 


± 


0.06 


1, 


,06 


± 


0.08 


1, 


,09 


± 


0.09 


1.05 


± 


0.07 


1, 


,09 


± 


0.10 


old 


NGC 


6341 


-2, 


,02 


-1, 


.94 


3.72 


± 


0.04 


1, 


.02 


± 


0.04 


1, 


.03 


± 


0.04 


1, 


,03 


± 


0.05 


1.05 


± 


0.04 


1, 


.03 


± 


0.05 


old 


NGC 


6352 


-0, 


,36 


-0 


.56 


4.13 


± 


0.07 





.90 


± 


0.06 





.99 


± 


0.07 


1, 


,03 


± 


0.08 


0.99 


± 


0.07 


1, 


.02 


± 


0.10 


old 


NGC 


6362 


-0. 


,96 


-0, 


.85 


4.09 


± 


0.05 


1, 


.06 


± 


0.05 


1, 


.06 


± 


0.05 


1, 


,08 


± 


0.05 


1.04 


± 


0.03 


1, 


.07 


± 


0.06 


old 


NGC 


6366 


-0, 


,44 


-0 


.59 


4.16 


± 


0.12 





.95 


± 


0.07 


1, 


.04 


± 


0.13 


1, 


,08 


± 


0.16 


1.04 


± 


0.13 


1, 


.07 


± 


0.18 


old 


NGC 


6388 


-0, 


,60 


-0, 


.63 


4.05 


± 


0.08 


0, 


.90 


± 


0.07 


0, 


,94 




0.08 


0, 


,95 




0.10 


0.92 


± 


0.07 


0, 


.94 


± 


0.10 


old 


NGC 


6397 


-1, 


,72 


-1, 


.54 


3.79 


± 


0.04 


1, 


.01 


± 


0.04 


0, 


.99 


± 


0.04 


0, 


,99 


± 


0.04 


0.99 


± 


0.04 


0, 


.99 


± 


0.04 


old 



Table 4 — Continued 



Name [M/H]zw [M/H]cG M^g^ D07zw D07cg P04cg B94cg GOOcg Group 



NGC 6441 


-0.45 


-0.46 


4.04 


± 


0.08 


0.85 


± 


0.06 


0.88 




0.07 


0.90 




0.07 


0.88 




0.05 


0.85 


± 


0.09 


old 


NGC 6496 


-0.36 


-0.56 


4.11 


± 


0.07 


0.88 


± 


0.06 


0.97 


± 


0.07 


1.00 


± 


0.09 


0.95 


± 


0.07 


0.98 


± 


0.10 


old 


NGC 6535 


-1.56 


-1.29 


3.68 


± 


0.11 


0.87 


± 


0.04 


0.82 


± 


0.09 


0.82 


± 


0.09 


0.81 


± 


0.09 


0.83 


± 


0.08 


young 


NGC 6541 


-1.57 


-1.31 


3.88 


± 


0.04 


1.06 


± 


0.04 


1.01 


± 


0.04 


1.02 


± 


0.04 


1.00 


± 


0.03 


1.01 




0.04 


old 


NGC 6584 


-1.32 


-1.10 


3.82 


± 


0.03 


0.92 


± 


0.03 


0.88 


± 


0.03 


0.87 


± 


0.03 


0.89 


± 


0.03 


0.87 




0.03 


young 


NGC 6624 


-0.36 


-0.56 


4.12 


± 


0.07 


0.89 


± 


0.06 


0.98 


± 


0.07 


1.01 


± 


0.08 


0.98 


± 


0.07 


1.00 


± 


0.10 


old 


NGC 6637 


-0.58 


-0.64 


4.13 


± 


0.07 


0.96 


± 


0.07 


1.02 


± 


0.07 


1.07 


± 


0.08 


1.02 


± 


0.07 


1.05 


± 


0.09 


old 


NGC 6652 


-0.75 


-0.76 


4.07 


± 


0.06 


0.97 


± 


0.06 


1.01 


± 


0.06 


1.03 




0.08 


1.01 


± 


0.06 


1.03 


± 


0.08 


old 


NGC 6656 


-1.53 


-1.27 


3.87 


± 


0.05 


1.03 


± 


0.05 


0.99 


± 


0.05 


0.99 


± 


0.05 


0.98 


± 


0.05 


0.98 


± 


0.05 


old 


NGC 6681 


-1.42 


-1.13 


3.93 


± 


0.04 


1.07 


± 


0.04 


1.00 


± 


0.04 


1.01 


± 


0.04 


1.00 


± 


0.03 


1.01 


± 


0.04 


old 


NGC 6715 


-1.32 


-1.03 


3.79 


± 


0.03 


0.90 


± 


0.03 


0.84 


± 


0.03 


0.83 


± 


0.03 


0.83 


± 


0.04 


0.84 


± 


0.02 


young 


NGC 6717 


-1.11 


-0.87 


4.06 


± 


0.05 


1.09 


± 


0.05 


1.03 


± 


0.05 


1.07 


± 


0.05 


1.03 


± 


0.04 


1.04 


± 


0.06 


old 


NGC 6723 


-0.90 


-0.82 


4.07 


± 


0.05 


1.02 


± 


0.05 


1.02 


± 


0.05 


1.06 


± 


0.05 


1.02 


± 


0.04 


1.04 


± 


0.06 


old 


NGC 6752 


-1.32 


-1.02 


3.88 


± 


0.03 


0.98 


± 


0.03 


0.92 


± 


0.04 


0.91 


± 


0.03 


0.91 


± 


0.02 


0.92 


± 


0.04 


young 


NGC 6779 


-1.72 


-1.50 


3.89 


± 


0.05 


1.10 


± 


0.05 


1.07 


± 


0.05 


1.06 


± 


0.06 


1.05 


± 


0.04 


1.06 


± 


0.04 


old 


NGC 6809 


-1.58 


-1.32 


3.83 


± 


0.04 


1.01 


± 


0.04 


0.96 


± 


0.04 


0.96 


± 


0.05 


0.95 


± 


0.04 


0.96 


± 


0.05 


old 


NGC 6838 


-0.44 


-0.59 


4.19 


± 


0.07 


0.98 


± 


0.08 


1.07 


± 


0.08 


1.10 


± 


0.10 


1.06 


± 


0.08 


1.11 


± 


0.10 


old 


NGC 6934 


-1.32 


-1.10 


3.81 


± 


0.04 


0.91 


± 


0.04 


0.87 


± 


0.04 


0.86 


± 


0.04 


0.88 


± 


0.04 


0.86 


± 


0.03 


young 


NGC 6981 


-1.28 


-0.99 


3.82 


± 


0.03 


0.91 


± 


0.02 


0.85 


± 


0.02 


0.84 




0.02 


0.85 




0.03 


0.85 


± 


0.02 


young 


NGC 7078 


-1.91 


-1.80 


3.74 


± 


0.05 


1.00 


± 


0.04 


1.01 


± 


0.04 


1.01 




0.04 


1.02 




0.04 


1.01 


± 


0.04 


old 


NGC 7089 


-1.39 


-1.09 


3.86 


± 


0.04 


0.98 


± 


0.04 


0.92 


± 


0.04 


0.91 


± 


0.04 


0.91 


± 


0.03 


0.92 


± 


0.05 


young 


NGC 7099 


-1.83 


-1.70 


3.76 


± 


0.05 


1.01 


± 


0.04 


1.01 


± 


0.05 


1.01 


± 


0.04 


1.01 
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